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Dyson processes on the octonion algebra 


Songzi Li 1 


Abstract We consider Brownian motion on symmetric matrices of octonions, 
and study the law of the spectrum. Due to the fact that the octonion algebra is 
nonassociative, the dimension of the matrices plays a special role. We provide 
two specific models on octonions, which give some indication of the relation 
between the multiplicity of eigenvalues and the exponent in the law of the 
spectrum. 


1 Introduction 

The study of the laws of the spectrum is one of the most important topics 
in random matrix theory. One may consider stochastic diffusion processes 
on specific set of matrices, for example symmetric or Hermitian matrices. 
Usually one considers the empirical measure of the spectrum, which is often 
again a stochastic diffusion process, called a Dyson process, see the works of 
Wigner [IS], Mehta na. Dyson [7j, Anderson-Guionnet-Zeitouni [I], Erdos 
and co 011311 ] , Forrester m and references therein. 

Let us recall some classical results on this topic. Consider specific matri¬ 
ces with independent Gaussian elements : real symmetric (/3 = 1), Hermi¬ 
tian (/3 = 2) and real quaternionic (/3 = 4). Then the law of their eigenval¬ 
ues (Ai)i<j< n , ordered as Ai < ... < A„, has a density with respect to the 
Lebesgue measure dA\...dA n which is 



( 1 ) 


j =1 l<j</c<n 
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where Cp n are constants depending on 0, n. 

On the other hand, one can also consider random matrices with stochastic 
process as entries. In his paper [7], Dyson derived the stochastic equations 
of the eigenvalues of Hermitian matrices whose elements are independent 
complex Brownian motions (see also in this direction Anderson-Guionnet- 
Zeitouni [T], IVIeht a 111 i , Li-Li-Xie 131). The stochastic process on the spec¬ 
trum provides a dynamic way to study the law of the eigenvalues of the 
matrix with Gaussian entries. In fact there are two ways: one is through the 
law of the eigenvalues of matrices with Brownian motions as entries, con¬ 
sidered at time t = 1; the other one is through the matrix whose elements 
are Ornstein-Uhlenbeck process, since when t —> oo the law of the matrix 
converges to a matrix with Gaussian entries, and the law of its spectrum is 
invariant through O-U process: in this case, the law of the spectrum may 
be seen as the invariant (in fact reversible, see definition U) law of the pro¬ 
cess. This is in general a much easier way to identify the law, since reversible 
measures are easy to identify through the knowledge of the generator. 

Meanwhile, if we consider real symmetric matrix, Hermitian matrix and 
real quaternionics matrix as real ones of dimension respectively nxn, 2nx 2 n, 
4 n x 4n, the multiplicity of the eigenvalues is again 1, 2 and 4. This fact 
leads us to wonder whether this exponent factor in the density reflects the 
multiplicity of the eigenvalues. However, this is not true. 

In a recent paper, Bakry and Zani [4], the authors considered real symmet¬ 
ric matrices whose elements are independent Brownian motions depending on 
some associative algebra structure of the Clifford type. Their computation 
of the law of the spectrum shows that, even though there is still the term 
ni< i<fc<n |A& — Aj|^ with /3 = 1,2,4, the factor 0 here reflects the structure 
of the algebra, known as Bott periodicity, rather than the dimension of the 
eigenspaces, which in this situation may be as large as we want. 

The previous study on Dyson Brownian motion, including the work of 
Bakry and Zani [3] on Clifford algebra, mainly concentrated on the case where 
the underlying algebra is associative. It is therefore worth understanding how 
important this property is in the study of the related Dyson processes. The 
octonion algebra, which is nonassociative but only alternative, provides a 
good example for us to start with. Its structure differs from the Clifford 
one, although Clifford algebras with 1 or 2 generators coincide with complex 
numbers and quaternions, the Clifford algebra with 3 generators does not 
coincide with octonions, even if the algebras have the same real dimension 8. 
In his book nn, section 1.3.5, Forrester mentions that the distribution Q 
with 0 = 8 can be realized by 2 x 2 matrices on octonions, with Gaussian 
entries. It is therefore worth to look at the associated Dyson process, which 
could also provide this result through the study of its reversible measure. 

There are only four normed division algebras: R, C, H and O. We are 
familiar with R, C, and while the quaternion algebra H is noncommutative 
but associative,the octonion algebra O is nonassociative, but only alternative. 
Even though their properties are not so nice, octonions have some important 
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connections to different fields of mathematics, such as geometry, topology and 
algebra. One interesting example is its role in the classification of simple Lie 
algebra. There are 3 infinite families of simple Lie algebras, coming from the 
isometry groups of the projective spaces RP", CP" and HP”. The remaining 5 
simple Lie algebras were later discovered to be in connection with octonions: 
they come from the isometry groups of the projective planes over O, O <8> C, 
0®H, © Cg> O and the automorphism group of octonions. It is also worth to 
mention that, according to the independent work by Kervaire jT2j and Bott- 
Milnor [Q in 1958, there are only 4 parallelizable spheres: §°, S 1 , § 3 and § 7 , 
which correspond precisely to elements of unit norm in the normed division 
algebras of the real numbers, complex numbers, quaternions, and octonions. 
See more examples in the paper by J. Baez [2j. 

For the eigenvalue problem of matrices on octonions, Y.G.Tian proved in 
his paper m that 2x2 Hermitian matrix on octonions has 2 eigenvalues, 
each of them has multiplicity 8. For 3x3 Hermitian octonionic matrix, Dray- 
Manogue [6] and Okubo m showed that it has 6 eigenvalues with multiplicity 
4. For 4x4 and 5x5 Hermitian octonionic matrices, there are only numerical 
results, indicating that the eigenvalues have multiplicity 2 m- It is still 
unknown for matrices in higher dimension. Following the analysis of Bakry 
and Zani U, one may expect that the study of probabilistic models on 
matrices of octonions could give new insights in these directions. 

In this paper, we consider Brownian motions on symmetric matrices of 
octonions. Due to the fact that octonions are nonassociative, and in contrast 
with the Clifford case, the dimension of the matrices plays a specific role. In 
fact, contrary to the real, complex and quaternionic cases, octonions do not 
give rise to infinite series of Lie groups but only specific ones, which are closely 
related to dimension 2. Thus the study of Dyson processes is mainly pertinent 
in this dimension, although we introduce another probabilistic model related 
to the octonion algebra, but with a special structure, see Section 4.2. To 
study the law of the spectrum of the matrices, we consider the processes on 
the characteristic polynomials P(X), as introduced in the paper by Bakry 
and Zani [4]. Because of the specific structure of octonions, the traditional 
way to compute the law of the spectrum turns out to be quite hard, while 
computation on the process of P(X) provides a simpler and more efficient 
method to see things clearly. 

The paper is organized as follows. Section 2 gives an introduction to the 
basics of the octonion algebra; Section 3 explains briefly the language and 
tools of symmetric diffusion process; Section 4 states our main results, two 
specific models on octonions, and then we explain what is so special about 
dimension 2; Section 5 is devoted to the demonstration of the connection 
between the algebra structure and the Euclidean structure associated with 
the associated symmetric matrices, and the fact that the the two exponents, 
multiplicity of eigenvalues and exponent in the law of eigenvalues, are not 
correlated. 


4 


Songzi Li 


2 The octonion algebra 

In this section, we recall some facts about the octonion algebra, and we refer 
to (2] for more details. We start with a few definitions. 

Definition 1. An algebra A is a division algebra if for any a,b £ A, with 
ab = 0, then either a = 0 or b = 0. A norrned division algebra is a division 
algebra that is also a norrned vector space with ||a6|| = ||a||||6||. 

Definition 2. An algebra A is alternative if the subalgebra generated by any 
two elements is associative. By a theorem of Artin m, this is equivalent to 
the fact that for any a, & € A, ( aa)b = a(ab), (ba)a = b(aa). 

As mentioned earlier, there are only four norrned division algebras, R, C, 
H, O. There is a nice way called "Cayley-Dickson construction" to produce 
this sequence of algebras: the complex number a + ib can be seen as a pair 
of real numbers (a, &); the quaternions can be defined as a pair of complex 
number; and similarly the octonions is a pair of quaternions. As the con¬ 
struction proceeds, the property of the algebra becomes worse and worse: 
the quaternions are noncommutative but associative, while the octonions are 
only alternative but not associative. 

Since octonions and Clifford algebra are both the algebra with dimension 
2™ (in this case n = 3), which share some special property, we can use the pre¬ 
sentation provided in Bakry and Zani [4] to describe the algebra structure on 
a basis of octonions, in order to simplify the computations. This presentation 
is not classical, and we shall therefore use the table below. 

Define E = (1, 2, 3}, and let V(E) denote the set of the subsets of E. For 
every set A £ V(E),we associate a basis element u>a in the octonion algebra, 
with = Id, the identity element. Then an element x £ O can be written 
in the form 

x = E a 'a^a, xa £ R, 

A 

and the product of two elements x and y is given by 

xy = '^2 xaVbuaub- 

a,b 

It remains to define loa^b for given A,B£ V(E), which is given through 
the following rule: denote by A.B the symmetric difference A U B \ {A n B), 
then u>aub = ( A\B)uja.b , where (A\B) takes value in { — 1,1}. Then, the 
multiplication rule in the octonion algebra is defined by a sign table, which 
is as follows : 
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0 

{1} 

{2} 

{3} 

{1,2} 

{1,3} 

{2,3} 

{1,2,3} 

0 

1 

1 

1 

1 

1 

1 

1 

1 

{1} 

1 

-1 

1 

1 

-1 

-1 

1 

-1 

{2} 

1 

-1 

-1 

1 

1 

-1 

-1 

1 

{3} 

1 

-1 

-1 

-1 

1 

1 

1 

-1 

{1,2} 

1 

1 

-1 

-1 

-1 

-1 

1 

1 

{1,3} 

1 

1 

1 

-1 

1 

-1 

-1 

-1 

{2,3} 

1 

-1 

1 

-1 

-1 

1 

-1 

1 

{1,2,3} 

1 

1 

-1 

1 

-1 

1 

-1 

-1 


In this table, the element (i,j) is the sign (Ai\Aj), where A, is the ith element 
in the first column, Aj is the jth element in the first row. 

From the facts that for A, B ^ 0, = —1 and uja^b = —wbwa, it is 

easy to get the following rules: 


(A\A) 

(A\B) 


f - 1 ) A ± 0 , 

11, A = 

-(B\A), for B ^ A, A, B ^ 


It can be seen from the above table that O is an algebra, non-associative 
but alternative. Moreover O can be equipped with the Euclidean structure 
obtained by identifying O as a 8 dimensional (real) vector space via 

x = X A^A H► (^0, £{l}, £{2}, ^{3}, £{1,2}, £{1,3}, ^{2,3}, ^{1,2,3}) , 

A 

so that the inner product and the norm are respectively : 

{x,y) = J2 xa va , \\ x \\ = C^2 x a) 1/2 , 

A A 

so that {oja, A € V(E)} form a real orthonormal basis for the algebra O. 

Let us recall that to prove that O is a division algebra, it is usual to 
introduce the conjugate 

x = ^2 x a^a '-t x * = ^2 x aUa(A\A), 

A A 

and observe that (xy)* = y*x* , xx* = x*x and ||a;|| 2 = xx* , so that \\xy\\ 2 = 
(a :y){xy)* = (xy)(y*x*) = x{yy*)x* = ||x|| 2 ||j/|| 2 . 

Altough the previous table does not provide an associative algebra, the 
octonion algebra satisfies however some useful identities. In what follows, we 
shall make a strong use of Moufang identities, which are stated as follows : 
for elements x, y,z belongs to O, we have 
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z(x(zy)) = (zxz)y, 

(( xz)y)z = x(zyz), 
(zx)(yz) = (z(xy))z, 
(zx)(yz) = z((xy)z). 


We shall mainly use this for the elements uja G O, although Moufang 
identities provide more information than this. 

According to the alternativity property, we can get some basic formulae 
about the sign table {(A|£?)} for octonions, 

Lemma 1. For A, B,C,D G V(E), we have 

1. ( A.B\B) = (A\B)(B\B). 

2. (A.B\A)(A.B\B) = (A.B\A.B). 

3. If A.B ± 0, (A.C\A){B.C\B) = ~(A.C\B)(B.C\A). 

4. If A.B.C.D = 0, 

{B.C\C)(C.D\D)(D.A\A)(A.B\B) = (B.D\B.D). 

Proof. The first one is just the result of alternativity: 


while 

Hence 


(loaixb)ixb = ( A\B)uja.b<xb = (• A.B\B)u a , 


(ujau)b)ub = uj a (u)b<xb) = ( B\B)w > a , 


(A.B\B) = (A\B)(B\B). 


The second one can be easily proved by the first statement. 

For the third statement, we first remark that, for any vector x = xc^Ci 
xu }a is always orthogonal to x ujb if A ^ B. Indeed, to see this, we may reduce 
to the case where ||x|| = 1, and then observe that the fact that the algebra 
is a division algebra shows that for any y £ O, y xy is an orthogonal 
transformation. For A,Bg V{E ), A.B ^ 0, and any C G V{E), choose 
D = A.B.C and set x = u>c + wd- Then 


0 = ( XU]A, XU}B) 

= (( Cj A)uic.a + {D\A)ujd.a, ( C\B)ujc.b + (D\B)ujd.b) 
= {{C\A)ujc.a, ( D\B)ujd.b) + {{D\A)wd.a, (C\B)ujc.b) 
= (C\A)(D\B) + (D\A)(C\B). 


Since D = A.B.C, the above formula indicates that for A,B,C G 'P(E) 
A.B ± 0, 
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(C\A)(A.B.C\B) + (A.B.C\A)(C\B) = 0. 

By changing C into AC, we get 

(A.C\A)(B.C\B) + (B.C\A)(A.C\B) = 0, 

then the third statement is proved. 

For the last statement, denote 

0 := (B.C\C)(C.D\D)(D.A\A)(A.B\B). (2) 

Then, from Moufang identities, 

((wb.cwc)(whjwa))((wc.iiwd)(w j 4.bwb)) = (B.C\C)(C.D\D)(D.A\A)(A.B\B)(uibud){ucva) 

= 0(B\D)(C\A)lob.d^Jc.a 
= 0(B\D)(C\A)u>b.dub.d 
= G{B\D)(C\A)(B.D\B.D). 


On the other hand, 

ojaojb = {A.B\A.B)(A\A)(B\B)cjbwa- 


{{wb.CUc)(ud.AUa))((uc.DUd)(ua.BWb)) 

= (D\D)(B.C\B.C)(A\A)(A\A)(B\B)(A.B\A.B)((ujb.cuc){uaub.c))((ua.bud){ubua.b)) 
= ( D\D)(B.C\B.C)(B\B)(A.B\A.B)(ujb.c{ucua)ub.c){ua.b(udub)ua.b) 

= (D\D)(B.C\B.C)(B\B)(A.B\A.B)(C\A)(D\B)(u J b.cu>c.aUb.c)(ua.bUd.bUa.b) 

= (D\D){B\B)(C\A)(D\B)w c .awd.b 
= (D\D)(B\B)(C\A)(D\B)(B.D\B.D). 


Hence, 

0 = (B\D)(C\A){B.D\B.D){D\D)(B\B){C\A)(D\B)(B.D\B.D) 
= (B\D)(D\D)(B\B)(D\B) 

= (B.D\B.D), 


which ends the proof of the lemma. 

For a nxn matrix on octonions, write it as M. = ^4 M a u)a, where {M' 4 } 
are real nxn matrices. For an n dimensional vector X b ojb, 
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C^ j M a u) A )C^ j X b uj b ) = ^M a X b (A\B)<jo a .b = J2( A -B\B)M a - b X b uj a . 

A A,B A,B 

Therefore, At can be expressed by the real 8 n x 8n block matrix {M A ' B }, 
where M A ’ B = (A.B\B)M ab . 

This leads to the following definition: 

Definition 3. A (2 3 x n) x (2 3 x n) block matrix M A ’ B (where A,B C 
{1,2,3}) is a real octonionic if M A ’ B = (A.B\B)M ab , where M A = M A ’® 
is a family of 8 n x n square matrices. It is the real form of a matrix with 
octonionic entries. We shall denote it as At = Y1,a M a loa- 

Then, we shall say that an octonionic matrix is symmetric if its real form 
is symmetric. This corresponds to the fact that, for any A £ V(E), ( M A y = 

(. A\A)M a . 

That is to say, (M A ’ B y = ( A.B\B)(M AB y = M B ’ A = ( B.A\A)M ab . 
Due to property 2 of Lemma 2.3, this leads to the fact that for any A £ V(E), 

( M A y = (A\A)M a 1 i.e. M® is symmetric while M A is antisymmetric for any 

A ^ 0. 

It is worth to point out that since the octonion algebra is not associative, 
there is no matrix presentation of the algebra structure for the octonions, and 
therefore the matrix multiplication of the real octonionic matrices does not 
corresponds to the octonionic multiplication of the associated matrices with 
octonion entries. Even the product of octonionic matrices is not octonionic 
in general. 

The inverse of an octonionic matrix is in general not octonionic, and its 
exact structure is not easy to decipher; the octonionic property may not be 
preserved. The following lemma gives a condition for this last property to 
hold, and will play an important role in the rest of this paper. 

Lemma 2. Let M = M a lo a be an octonionic matrix such that M® is 
invertible. Assume moreover that, for any A, B £ V(E) 

M a (M®)~ 1 M b = M b (M®)- 1 M a , (3) 

and that M c is invertible. Then, M is invertible and its in¬ 
verse N is octonionic, satisfying N = ^2 ^ujaX a , with 


N a = —7V 0 M a (M®) -1 , 

for A ^ 0, 

(4) 

IV® = M C (M®)“ 


(5) 


c 


Proof. In fact, assume the octonionic matrix N A,B = ( A.B\B)N ab = 
( A.B\A)N a b is the inverse of M , where 


N ab = ( A.B\A.B)N ab 


-N a b , A.B ^ 0; 
TV®, A.B = 0. 
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Then 

^2(A.C\A)(C.B\B)N a - c M cb = 0, for A ^ B, (6) 

c 

J2(A.C\A)(C.A\A)N a - c M ca = Id. (7) 

c 

Changing C into A.B.C in ®, we get Y,c( B - C \A)(A.C\B)N b - c M c - A = 0, 
then it is enough to have 

(A.C\A)(C.B\B)N a - c M cb + (B.C\A)(A.C\B)N b - c m c - a = 0. 


According to Lemma [lj it holds as soon as 

N ac M cb = n bc m c - a . 


( 8 ) 


Choosing C = B and then setting D = A.C , this leads to 

N d = N®M d (M®)-\ (9) 

for every D £ V(E). Now choose C = 0 in J5]) and apply (0 to N A and N B , 
plugging into 0, this reduces to 

7V 0 M a (M 0 )- 1 M b = AT® M b (M 0 )- 1 M a . 


Now 0 is 

N d M d = Id, 

D 

and using 0, this gives 

Af 0 (^M c (M 0 )- 1 M c )- 1 =Id, 

which means that IV 0 is invertible, and gives its inverse, such that © holds 
true. Then we can use © and 0 to get 0. 

Remark I. It is worth to observe for later use that if the matrix A4 on the 
octonions satisfies the assumptions of Lemma [21 then it is also the case of 
M - Xld. 


3 Symmetric diffusion operators on matrices 


We introduce the basics on symmetric diffusion operators, in a simplified 
version adapted to our case. For further details see [3j. 

Let E be an open set in R”, endowed with a er-finite measure n and let 
Ao be the set of smooth compactly supported functions, or of polynomials 
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functions on E. For any linear operator L : Aq i —> Ao, we define its carre du 
champ operator as 

r(f,g) = ±(L(fg)-fL(g)-gL(f)). 

We have the following 

Definition 4. A symmetric diffusion operator is a linear operator L: Ao®l K > 
Aq, such that 

1- L(l) = 0, 

2. \/f,g <E A 0 ®1, f fL(g) dfi = f gL(f) dg, 

3. v/ e A 0 ,r(f,f) > o, 

4 . V / = (/ i , •'' > fn), where /,; € Aq , <P is a smooth function R n i —> R and 
*(0) = 0, 


l m)) = E d&ifWi )+E /#)■ ( l0 ) 

* i,j 

Consider an open set El C E, and a given system of coordinates ( x *), then 
we can write 

L (/) = E 9 l3 {x)d^f + E h \x)dif, 

ij i 

where 

g^(x) = r(x l ,xi), b l {x) = L(P). 

In this paper, we perform computations on the characteristic polynomial 
P(X) = det (M — Xld) of a matrix M. Assume that we have some diffusion 
operator acting on the entries of a matrix M , described by the values of 
L (Mij) and M^i) for any ( i,j,k,l). Then, we have, 

P(logP(A),logP(r)) = E 5 Mw log(P(X))0M fe! log(P(y))P(M ij ,M w ), 

L(log P(X)) = E d M i5 log(P(X))L(M ij ) + E d MiJ d Mkl log(P(X))P(M ij , M kl ). 

i,j i,j,k,l 

To compute log (P(X)) and dM^M^ log(P(A)) in the above formulae, 

we use the Lemma 6.1 in Bakry and Zani [3], which we quote here without 
proof. 

Lemma 3. Let M = ( Mij) be a matrix and M~ l be its inverse, on the set 
{det-M ^ 0} we have 


d Mij log det M = M^ 1 , 
d Mij d Mkl log det M = 
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Hence, with M~ 1 {X) = (M - A/d)” 1 , 

P(logP(A),logP(Y)) = J2 M- 1 (X) ji M~ 1 (Y) lk r(M ij ,M kl ), (11) 

i,j,k,l 

L(logP(X)) = "^2 Mjl 1 (X)L(Mij) — J2 Mj k 1 (X)M l T 1 (X)r(M ij ,Mm 

i,j i,j,k,l 

According to Bakry and Zani [4] , one can get from P(P(A),P(Y)) and 
L(P(A)) informations about the multiplicities of the eigenvalues, and on the 
invariant measure of the operator L acting on P(A): 

If for some constants on, « 2 , cr 3 , 

L(P) = a±P" + 012 —p -, P(log P(X) ,log P(Y)) = ^^(^1 - ®). 

(13) 

and if there exists for some a £ M, a ^ 0 which satisfies 

a 2 (a i + a 2 ) — a(a i + 03 ) + a 3 = 0, (14) 


Then 

1. If a is a positive integer, it is the multiplicity of the eigenvalues of M; 

2. Write P(A) = n?=i(* — Xi) a , the invariant measure for the operator L 
in the Weyl chamber {x\ < ... < x n } is, up to a multiplicative constant, 

dn = ([[{xi-Xj) ) “3 duo, 

i<j 

where d/i 0 is the Lebesgue measure. 


4 Symmetric matrices on octonions 

Our aim is to describe the law of the spectrum of the real form of symmetric 
matrices on octonions. The block matrix is At = ((A.B\B)M ab )a t B^v(E)i 
satisfying ( M A ) 1 = ( A\A)M a from the symmetry assumption. 

We will focus on cases where the symmetry condition d3j) of matrix At — 
Aid is satisfied, i.e. where the matrix 

U{X) ■= (At - Aid )" 1 

is octonionic (almost surely for the stochastic process under consideration). 
Setting 

P(A) := det(At - Aid), 


by Lemma [3] we have 
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r(\ogP(X),\ogP{Y)) 


E 


u ?i’ Au ik’ c r(M?j 


A,B 


A,B,C,D 

i,j,k,l 




) 


E {A.B\A.B){C.D\C.D)U£ B TJg D r{M$ B , M$S) 

A,B,C,D 

i,j,k,l 


(where we used property 2 of Lemma 2.3), and 


L(logP(X))=Y,U A ’ B (X)L(M£ B )- E Uff{X)uS' A 

A,B A,B,C,D 

A A A A b* 1 

(x)r(M£ B ,M' t 

(16) 

For further use in both examples, we state (without proof) two preliminary 
lemmas. The first lemma collects some elementary facts, consequences of the 
definition, the symmetries and property 2 of Lemma |T| 

Lemma 4. For F £ V{E), we have 


trace U(X) F = Y / U( x )u, 

(17) 

trace U ( X) = 8 trace U (X) 0 , 

(18) 

trace [U(X) F U(Y) F ] = £ U(X) F U{Y) F , 

ij 

(19) 

{F\F)trace [U(X) F U{Y) F ] = £ U(X) F U{Y) F , 

(20) 


ij 


trace \U{X)U{Y)] = 8 ^ ^{C\C)trace [U(X) C U{Y) C ] . (21) 

c 

The second lemma gives expressions of the traces in terms of the characteristic 
polynomials. The first two identities are obtained by derivation from 


logP(X) = trace \og{M — Xld) = —trace log[/(X), 
and the third one is a consequence of the first one and the resolvent equation. 

Lemma 5. 


trace U(X) 


trace ( U(X ) 2 ) 
trace [U{X)U(Y)} 


pvn 

P(X) ’ 

P'(X) 2 P"(X) 

P(X) 2 P(X) ’ 

1 fP'(X) P'(Y)\ 
Y-X \P(X) P(Y) ) ' 


( 22 ) 

(23) 

(24) 
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4-1 The dimension 2 case 


Consider M = M a lja, where {M A } are matrices whose elements are in¬ 
dependent Brownian motions. For A ^ 0, due to symmetry of A4, (M A ) 4 = 
( A\A)M a = —M a . Such matrices naturally satisfy the symmetry restric¬ 
tion [3] in dimension 2, since the 2x2 antisymmetric matrices are all of the 

form 0~)’ an< ^ ^ e y are therefore all proportional to each other. This 

property fails in higher dimensions. Set 

M f ) = \sA, B {SikSii + {A\A)S a S jk ), L(M a ) = 0, (25) 


which reflects the symmetry of the matrices. Notice that the inverse matrix 
U{X) is also symmetric with (U A ) t = (A\A)U a . We have the following result 


Proposition 1. For the 2x2 symmetric matrix M = ^ M a uja, 


r(\ogP(X),\ogP(Y)) = (^1 


P'(Y) \ 
P(Y) ) ’ 


P(logP(X)) = 3 


/P'(X) 2 

Iw 


P"{ X) \ 
P{X) ) 


1 P'(X) 2 

2 p(xy ' 


(26) 


(27) 


Proof. Proof of (l26l) : From (fl5l) and (l25|) we have 


P(logP(W), log P(Y)) 

= E {A.B\A.B){C.D\C.D)U{X)%- A U{Y)^ c \{8 ik 8 jl + (A.B\A.B)5 u 6 jk ) 

A,B,C,D 

AB=CD 


= 8j2tr(U(X) F U(Y) F )(F\F) 


= 8tr (U(X)U(Y)) = 


8 fP'(X) 


Y-X 


\P(X) 


P'(Y) \ 
P(Y) ) ’ 


where we applied EOjl . (fTH . (ET1) and formula 3 in Lemma [5] This ends the 
proof of EH . 

Proof of E71) : In the following, we write U is for U{X). From (TIH) and 
(l25l) . we have 

L(logPpO) = -i E (B-D^.D^iUff 0 ) 2 

Z B,C,D ij 

E (B.D\B.D)(A.B\A.B)(J2u fc ) 2 . 

A,B,C,D i 

A.B.C.D^H) 


1 

2 
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On the one hand, in view of (EU1) and m 

E (B.D\B.D)J2 (U* C ) 2 = -6 trace ( U 2 ) . 

B,C,D ij 

On the other hand, since for F ^ 0, U F is antisymmetric, so that tr U F = 0, 

E (B.D\B.D)(A.B\A.B)(J2 U u C f = 8 2 (trace t/ 0 )*. 

A,B,C,D i 

A.B.C.D=<D 

Combined as in (128ft . these sums give : 

L(log P(X)) = 3tr U 2 - 4 • 8 (tr (7 0 ) 2 = 3tr U 2 - | (tr U) 2 , (28) 

which, in view of ( 1221 ) and (E^l) ends the proof of C71) and then the proof of 
the Proposition. 

Remark 2. s By the results of above proposition and formula (flOl) . it is easy 
to get 

1 P'( Yl 2 

Now chose a\ = —11, a 2 = 11 — 5 and = 8 in formula (fl3l) : the resulting 
value for a is a = 8 . This shows that the multiplicity of the eigenvalues is 8 . 
Assume p is the density of the invariant measure of L of the coordinates 
{Xi } in the Weyl chamber, then according to our discussion in the previous 
section, we have 

p = C JJ(xj - Xj) 8 . 
i<j 


4-2 Another model in any dimension 


We now provide another set of random octonionic matrices for which the 
symmetry condition ([3J is automatically satisfied. 

Let iW® be a symmetric matrix with independent Brownian motions as its 
entries. For all A, B ^ 0, let M A = M B = A be a random antisymmetric 
matrix with independent Brownian motion as its off diagonal entries. Then 
consider M = + A^2 c ^ijJc- This model is similar to the Hermitian 

case considered in Bakry and Zani [J] (see Remark 0]). Similarly to the Her¬ 
mitian case, we set 




( 29 ) 
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r(Aij,Aki) = Y^(S lk Sji - 6uS jk ), (30) 

r(M? j ,A kl ) = 0 , L(M$) = o. (31) 

Due to Lemma[TJ for the inverse matrix U(X) = (M — X I) -1 , we have for 
every C/0 

U c = -[/ 0 M C (M 0 - XI)- 1 = -U^AiM® - XI)- 1 , 
which means for all C ^ 0, U c is the same, and we denote it by U a - 
Proposition 2. For the matrix M = M^uiq + Auc on the octonions, 


P(logP(X),logP(Y)) 

LQogP) 


8 fP'(X) 
Y-X V P(X) 
1P'(X) 2 
~8 P(X ) 2 ' 


P'(Y) \ 

P(Y) ) ’ 


Proof. On the one hand, from m and m 

r(log P(X), log P(Y)) = isr + ^S 2 , (32) 

where 


Si~ E E U(X)%U(Y)l(5 ik 6 jl + 6 a 6 jk ), 

A=B,C=Di,j,k,l 

S2 '-= E E U(X)f i B U(Y)g D (S lk S jl - 6 a 6 jk ). 

AjtB,Cj£D 


A careful computation, using the fact that [7(F ) 0 is symmetric and U(Y) a 
is antisymmetric gives 


Si 


2 - 8 2 trace 


P(X) 0 P(F ) 0 


—2-7 2 -8 2 trace [U(X) a U(Y) a ] , (33) 


and then, using (15^1) 


P(logP(X),logP(F)) 


8 2 trace 


U(X)*U(yY 


—7-8 2 trace [U(X) a U(Y) a \. 


Going back to m we see that 


tr (U(X)U(Y)) 


8 trace 


U(X)*U(Yf 


7 • 8 trace [U(X) a U(Y) a )( 34) 


so that 
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F(logP(X),logP(F)) = 8 tr [U(X)U(Y)] = 


On the other hand, with U for U(X) 

L(logP) = ~ E ( B.C\C)(DA\A)(A.B\B)(C.D\D)U(X)f k c U(X)^ A (5 lk S jl + SuS jk ) 

Z A—B.C—D 

~ E { B -C\C){D.A\A){A.B\B){C.D\D)U{X) B k c U{X) BA {5 lk 5 ]l -8 a 5 3k ) 

4 A^B,C^D 
i,j,k,l 



Let us first remark that 


E Uf k U%(5 x k6ji ± SuSjk) = (P|P)trace (U F U G ) ± (trace Z7 F )(trace U G ). 

ijkl 

(35) 


For the first part, we have 

S[ = Etrace ( U BC ) 2 + Y ( B - C \ B - C )[trace U BC ] 2 

B,C B,C 


= 8 |trace (P 0 ) 2 + 7trace ( U a ) 2 + 8 [trace P 0 ] 2 . 

For the second part, going back to the notation @ for <9 and applying 
(1551) we have 


S' 2 = E ° x ((-B.C|B.C , )trace [p^-Cp^] - trace P BC trace ^>6) 

A^B,C^D 

Let us split the sum into four parts according to B = C or not, and A = D 
or not. 

i) When A ^ D, B = C, the sum vanishes. Indeed, in this case, 

(. B.C\B.C)trace [U BC U AD ] —{trace U B C trace U A D ) = trace [L/®C/ a ], (37) 

and <9 = {A.D\A){A.B\B){B,D\D) which is antisymmetric in A,B. 

ii) The same occurs when A = D and B ^ C. 

iii) When A = D and B = C, 

{B ,C\B ,C)tr ace [U B - c U AD }-{trace P BC trace U A ' D ) = trace [P 0 ] 2 -(trace P 0 ) 2 , 

and <9 = (B.D\B.D) in view of property 4 of Lemma 2.3, so that the contri¬ 
bution is 
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(B.D\B.D) ((B.C\B.C)tiace [U B C U A - D ] - trace U B C trace U A D ) 

AytB,CjtD,B=C,A=D 

= -7 ■ 8 (trace [E / 0 ] 2 - (trace E/ 0 ) 2 ) . 
iv) When A^D,B^C, 

(B.C\B.C)tiace [U B C U A ' D ] — (trace E/ BC trace U A D ) = —trace ( U a ) 2 ■ 

With the help of some computer algebra, we get 

Y O = 2 3 • 7 2 , 

A^B,C^D,B^C,A^D 

and finally, all the contributions in (15U1) give so that 

S' 2 = — 7 • 8 (trace [E / 0 ] 2 — (trace E/ 0 ) 2 ^ — 7 2 • 8 trace (E / a ) 2 • 

Going back to (l35l) and (l36l) . we conclude, using again (IdTl) 

L(logP) = — 8 (trace E / 0 ) 2 = — i(trace U ) 2 , 

which ends the proof of the proposition. 

Remark 3. Similarly we have 

L(P(X)) = (8 - ~ 8 P"(X). 

Chose = — 8 , «2 = 8 — «3 = 8 , we still obtain the multiplicity a = 8 , 

while the density of the invariant measure of L is then 

cni *- x ^ 2 ' 

i<j 

Remark 4- Recall that in Bakry and Zani 0], section 7.1, for a Hermitian 
matrix H — M+iA with independent Brownian motions as its entries (where 
M is symmetric, A is anti-symmetric), we have 

= -( S ik 6ji + SuSjk ), 

r(Aij,A u ) = -{S ik Sji $u$jk)- 

In our model M = M 0 o >0 + wp, denote e the specific element in 

the octonion algebra, e = Notice that 
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e 2 = —7, 

which indicates that e works like i in the Hermitian matrices, just with a 
different variance. Therefore, this example is indeed similar to the case of 
Hermitian matrices. 

Remark 5. In this remark we would like to discuss why the dimension 2 is so 
special. Consider a more general model: let 

M = M 0 w o + mCu} c , 

C #0 

with M c = xc-Aq, where M 0 is a Brownian motion on symmetric matri¬ 
ces, {xc} a series of Brownian motions on R, and Aq = {atj} a fixed anti¬ 
symmetric matrix. Obviously this model satisfies the symmetry condition [3j 
When M is a 2 x 2 matrix, it can be considered as a special example of the 
first case. Let e = X C^C- Different from the previous model, in this 

case e can be considered as a Brownian motion on the basis of octonions 
satisfying e 2 = - Ec #0 \ x c? = H e | 2 - 
Therefore, 

= -(SikSji + SuSjk ), 

r(M$, M B ) = S A =Ba ij a k i,A , P ^ 0, 
r(M®.,M^) = O,H^0. 


Similar computations yield 

P(logP(X),logP(Y)) =8 ^ tr(P(X) 0 t/(Y) 0 )+8 ^ (tr{U(X) BA A 0 )tr{U(Y) BA A 0 )), 

A.B =0 A.B^H 

P(logP) = ~ J2 tr (U(X) A C ) 2 - 4tr([/(X) 0 ) 2 - 4(trP(X) 0 ) 2 
Z A.C #0 

+5 ^2 tr{U(X) BC A 0 U{X) BC A 0 ) + 56tr{U(X) l!, AoU(X^A 0 ), 
B.C^Q 


which are hard to describe in terms of P. When the matrix is 2 x 2, the 
following equalities hold: for C ^ 0, 

tr(U(X) c AoMU(Y) c A 0 ) = -tr (U(X) C U(Y) C ), 
ti(U(X) c A 0 U(X) c A 0 ) = —tr {U{X) C )\ 
tr(t/(X) 0 H o t/(X)%) = tr ([/ 0 ) 2 - (tr[/ 0 ) 2 . 
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which give rise to the results in the first model. 

However, in higher dimensions, the above conditions are hard to satisfy. 

In fact when n = 2, it is enough to take Aq = ^ 0^' = z ‘ ^ 

formula 


C/ 0 = (M® + z 2 Ho(M 0 )-1Ho) _1 , 

U {1} = -2(7 0 H o (M 0 )-1. 


It is easily seen that isa2x2 antisymmetric matrix and can be written 
as 

U {1} = AA 0 , 

where A = z 2 _ d ^ M t) ■ Compare it with the expressions of [7® and U^, 
we have Aq = v ^_ z M 0 AoM 0 . Since A d = —this leads to (M 0 Ao) 2 = 
I = —det(M®)/, which is impossible to hold for any symmetric matrix 
M 0 in higher dimensions. This restriction insures the first two conditions, and 
the third one is proved by this and the fact that tr(M 2 ) — (trM ) 2 = — 2 det(M) 
holds in dimension 2 . 


5 Some remarks 


Our two models provide examples where the multiplicity of eigenvalues and 
the exponent fj in the law are not related, which is in accordance with the 
conclusion in Bakry and Zani (4], that the exponent reflects the structure 
of the algebra while the multiplicity of the eigenvalues is decided by the 
dimension of the eigenspaces. 

As we have seen, the octonionic structure of the matrix plays an important 
role. For higher dimension, the problem may be studied by our method if we 
know the structure of the inverse matrix, which is not necessarily octonionic. 
The main obstacle is still the non-associativity, which prevents any matrix 
presentation for octonionic multiplication. Let us recall that the 3x3 matrices 
on octonions have been studied by Dray and Manogue |5] and Okubo m 
using algebraic method, showing that there are 6 eigenvalues with multiplicity 
4. It it still an open problem to provide a probabilistic model in this case which 
would lead to this conclusion. 
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